Non-equilibrium quantum transport through a dissipative resonant level 



(N 

o 

(N 

> 
O 

in 



-a 

o 
o 



Chung-Hou Chung 1,2 , Karyn Le Hur 3 ' 4 , Gleb Finkelstein 5 , Matthias Vojta 6 , and Peter Woffle 7,8 
1 Department of Electrophysics, National Chiao-Tung University, HsinChu, Taiwan, R.O.C. 

2 Physics Division, National Center for Theoretical Sciences, HsinChu, Taiwan, R.O.C. 

3 Center for Theoretical Physics Ecole Polytechnique and CNRS 91128 Palaiseau, France 
4 Department of Physics and Applied Physics, Yale University, New Haven, CT, USA 

5 Department of Physics, Duke University, Durham, NC 27708, U.S.A. 
6 Institut fur Theoretische Physik, Technische Universitat Dresden, 01062 Dresden, Germany 
7 Institut fur Theorie der Kondensierten Materie, KIT, 76128 Karlsruhe, Germany 
8 Institut fur Nanotechnologie, KIT, 76021 Karlsruhe, Germany 
(Dated: November 19, 2012) 

The resonant-level model represents a paradigmatic quantum system which serves as a basis for 
many other quantum impurity models. We provide a comprehensive analysis of the non-equilibrium 
transport near a quantum phase transition in a spinless dissipative resonant-level model, extending 
earlier work [Phys. Rev. Lett. 102, 216803 (2009)]. A detailed derivation of a rigorous mapping 
of our system onto an effective Kondo model is presented. A controlled energy-dependent renor- 
malization group approach is applied to compute the non-equilibrium current in the presence of a 
finite bias voltage V . In the linear response regime V — > 0, the system exhibits as a function of 
the dissipative strength a localized-delocalized quantum transition of the Kosterlitz-Thouless (KT) 
type. We address fundamental issues of the non-equilibrium transport near the quantum phase 
transition: Does the bias voltage play the same role as temperature to smear out the transition? 
What is the scaling of the non-equilibrium conductance near the transition? At finite temperatures, 
we show that the conductance follows the equilibrium scaling for V < T, while it obeys a distinct 
non-equilibrium profile for V > T. We furthermore provide new signatures of the transition in the 
finite-frequency current noise and AC conductance via the recently developed Functional Renor- 
malization Group (FRG) approach. The generalization of our analysis to non-equilibrium transport 
through a resonant level coupled to two chiral Luttinger-liquid leads, generated by the fractional 
quantum Hall edge states, is discussed. Our work on dissipative resonant level has direct relevance 
to the experiments in a quantum dot coupled to resistive environment, such as H. Mebrahtu et at, 
Nature 488, 61, (2012). 

PACS numbers: 72.15.Qm,73.23.-b,03.65.Yz 
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I. INTRODUCTION 

Quantum phase transitions (QPTs)^^ which separate 
competing ground states represent generic phenomena in 
solid-state systems at zero temperature. The transition 
is frequently found to be continuous, often times giving 
rise to a quantum critical point. In the neighborhood of 
a quantum critical point of a metallic system the finite 
temperature properties as a rule show non-Fermi liquid 
behavior—. In recent years, quantum phase transitions at 
the nanoscale have attracted much attention 4- ! 5 . Much 
of the effort has been focused on the breakdown of the 
Kondo effect in transport of a quantum dot due to its cou- 
pling to a dissipative environment. However, relatively 
less is known about the corresponding out-of-equilibrium 
properties 1 ^ 23 ' 25 • . A finite bias voltage applied across 
a nanosystem is expected to smear out the equilibrium 
transition, but the current-induced decoherence might 
act quite differently as compared to thermal decoherence 
at finite temperature T, resulting in exotic behavior near 
the transition. 

Meanwhile, understanding the interplay of electron in- 
teractions and non-equilibrium effects in quantum sys- 
tems is one of the most challenging open questions in 
condensed matter physics. Many of the theoretical ap- 



proaches that have been proven so successful in treating 
strongly correlated systems in equilibrium are simply in- 
adequate once the system is out of equilibrium. The real- 
time Schwinger-Keldysh formalism 26 has been known 
as the most successful approach to non-equilibrium dy- 
namics since it offers a controlled perturbative expan- 
sion of the density operator. However, care must be 
taken to avoid the appearance of infrared divergences, 
in the perturbative approaches. Though much is known 
for quantum impurity systems in equilibrium, under- 
standing their properties in non-equilibrium steady-state 
is still limited. Nevertheless, significant progress has 
been made by different approaches, such as (1) analyti- 
cal approximations: perturbative renormalization group 
method (RG) 27 ' 28 , Hamiltonian flow equations 29 , Func- 
tional RG 30 ' 31 , strong-coupling expansions 3 ^, master 
equations^ 3 .; (2) exact analytical solutions: field theory 
techniques 3 ^, the scattering Bethe Ansatz^ 5 -, mapping of 
a steady-state non-equilibrium problem onto an effec- 
tive equilibrium system^ 6 -; (3) numerical methods: time- 
dependent density matrix renormalization group (RG) 37 , 
time-dependent numerical RG 38 , diagrammatic Monte 
Carlo^ 9 -, and imaginary-time nonequilibrium quantum 
Monte CarlrA 

In this paper, we provide a comprehensive analy- 
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sis of the non-equilibrium transport near a quantum 
phase transition in a dissipative resonant level model by 
employing the recently developed frequency-dependent 
RG2£ and Functional RG approaches^, and extending 
our earlier work in Ref. [lg. We aim to address several 
fundamental questions related to the non-equilibrium 
transport in quantum dot settings, such as: what is the 
distinct non-equilibrium conductance profile at zero tem- 
perature compared to that in equilibrium at finite tem- 
peratures near the transition? is there any scaling behav- 
ior of the conductance at finite temperatures and finite 
bias voltage near the transition? 

For this purpose, we investigate three classes of typical 
nano-models comprising a spinless resonant level coupled 
to: (i) two non-interacting Fermi-liquid leads subject to 
an Ohmic dissipative environment, where an Ohmic en- 
vironment can be realized in a nanoscale resistor and 
has many applications in phy sics ranging from meso- 
scopic physics (Refs. l8l lTfMll ) to biological systems^, 
(ii) two interacting fermion baths, in particular two Frac- 
tional Quantum Hall Edge (FQHE) 42 < leads, or the "chi- 
ral Luttinger liquids" where electrons on the edge of a 
2D fractional quantum Hall system show one-dimensional 
chiral Luttinger liquid behaviors with only one species 
of electrons (left or right movers), (iii) two interacting 
Luttinger-liquid leads subject to an Ohmic dissipative 
environment. 

In the class (i) model, the QPT separating the con- 
ducting and insulating phase for the level is solely driven 
by dissipation, which can be modeled by a bosonic bath. 
Dissipation-driven QPTs have been addressed theoret- 
ically and experimentally in various systems, such as: 
quantum dot systems*^, Josephson junction arrays^—, 
superconducting thin fil m 47 ' 48 , superconducting qubit^, 
qubits or resonant level systems coupled to photonic 
cavitie s 50 ' 51 , and biological system o 41 ' 52 . Here, we fo- 
cus on the non-equilibrium properties of the system near 
quantum phase transition. Meanwhile, for the class (ii) 
model, tunneling of electrons or quasi-particles between 
two FQHE states may in general suffer from the electron- 
electron interactions in FQHE. Interesting experimen- 
tally relevant questions arise regarding how interaction 
effects modify the noncquilibrium charge transport in 
such systems. Furthermore, one can extend the above 
two classes of models to a more general class (iii) model 
where both electron-electron interactions and the dis- 
sipation are present in the FQHE setpups, which have 
not been explored both theoretically and experimentally. 
Our results have relevance for the recent experiment in 
Ref . HH where the electronic transport through a resonant 
level in a nanotube exhibits the Luttinger liquid behav- 
ior, namely the conductance demonstrates a non-trivial 
power-law suppression as a function of bias voltage. 

This paper is organized as follows: In Section II A, 
the model Hamiltonian of class (i) is introduced. In Sec- 
tion II B, we establish rigorious mappings of our model 
system, the class (i) model at a finite bias voltage, onto 
the out-of-equilibrium anisotropic Kondo model as well 



as onto class (ii) and (iii) model systems subject to a 
finite voltage bias. We compute the current operator 
in Section II C for these three classes of models. We 
employ the nonequilibrium RG approach in Section III. 
Our results on nonequilibrium transport near the quan- 
tum phase transition both at zero and finite tempera- 
tures are presented in Section IV, followed by the results 
on the nonequilibrium finite-frequency current noise in 
Section V. We make a few remarks on the important is- 
sues of nonequilibrium quantum criticality in Section VI. 
Finally, we draw conclusions in Section VII. 

II. MODEL HAMILTONIAN 

A. Dissipative resonant level model 

Our Hamiltonian in all of the three classes of models 
mentioned above takes the following generic form: 

H = K^-^cliCki + Uc^d + h.c. (1) 

fc,i=l,2 

+ J2 X r(dU-l/2)(b r + bt) + Y,u r bib ri 

r r 

where U is the (real-valued) hopping amplitude between 
the lead i and the quantum dot, Cki and d are electron 
operators for the (Fermi-liquid type) leads and the quan- 
tum dot, respectively, /i^ = ±V/2 is the chemical poten- 
tial shift applied on the lead i (V will denote the bias 
voltage throughout this paper), while the dot level is at 
zero chemical potential. Here, b r are the boson opera- 
tors of the dissipative bath with an Ohmic type spectral 
density 8 -: J(w) = J2 r ^r^( w — = otu. Note that usu- 
ally we introduce a cutoff via a exp(— uj/u> c ) function in 
J {to)', here, we assume that w c is a large energy scale com- 
parable to the energy bandwidth of the reservoir leads. 
To simplify the discussion, we assume that the electron 
spins have been polarized through the application of a 
strong magnetic field. Note that our model can be real- 
ized experimentally in a quantum dot coupled to resistive 
environment as shown in Ref. |43| . 

In this section, we briefly summarize the behavior of 
our model system at equilibrium which means in the 
absence of a finite bias voltage (V — 0). A dissipa- 
tive resonant-level systems in equilibrium coupled to sev- 
eral leads maps onto the anisotropic one-channel Kondo 
model 8 - ' 10 ' 11 where the dimensionless transverse Kondo 
coupling gy is proportional to the hopping t between 
the level and the leads and the longitudinal coupling 
qf^ oc 1 — \J~ol (the exact prefactors are given in Refs. 

see also Sec. II B and Appendix A). Here, the 
superscript (e) in gfl z refers to the equilibrium cou- 
plings. The model exhibits a Kosterlitz-Thouless (KT) 
QPT from a delocalized (Kondo screened) phase for 
g^ + gf^ > 0, with a large conductance, G ~ e?/h, 
to a localized (local moment) phase for g 1 ^ + gf^ < 0, 
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with a small conductance, as the dissipation strength is 
increased (see Fig. 1). For g± — > 0, the KT transition 
occurs at a c = 1. As a —¥ a c , the Kondo temperature Tk 
obeysi lnT^f cx I /(a — a c ). Note that here we assume 
our resonant level system exhibits the particle-hole (p-h) 
symmetry; namely, the resonant-level energy ed is set to 
be zero (e^ = 0). However, in a more general resonant- 
level model where p-h symmetry is abscent, an additional 
term of the form e^d is present in the Hamiltonian Eq. 
@ • In terms of its equivalent Kondo model, this p-h sym- 
metry breaking term plays the role as an effective local 
magnetic field B z cx acting on the impurity spin in the 
Kondo modeUi, which needs more involved treatments 
and exceeds the scope of a simple and generic model sys- 
tem considered in the present work. 

In equilibrium, the dimensionless scaling functions 
gf^(T) and g z e \T) at the transition are obtained via the 
renormalization-group (RG) equations of the anisotropic 
Kondo model: 



^l(T) = - 5 il(T) = (21n(2VT)r 



(2) 



where V = Doe 1 ^ 29 ^, with Do being the ultraviolet 
cutoff. Having in mind a quantum dot at resonance, 
D = min(<5e, ui c ), with Se being the level spacing on 
the dot and w c the cut-off of the bosonic bath; Do is 
of the order of a few Kelvins. At low temperatures 
T <C -Do, the conductance drops abruptly with decreas- 
ing temperatures^: 



G eq (a c ,T < D ) cx 



9±jcr CO 



1 



ln z {T/V) 



(3) 



Below, we analyze the non-equilibrium (V ^ 0) 
transport of our model system at the KT transition and 
in the localized phase in the double-barrier resonant 
tunneling regime where the dissipative resonant level 
couples symmetrically to the two leads (ti = t 2 = t). 
Note, however, that when the dissipative resonant level 
couples asymmetrically to the leads t\ ^ ti, as has 
been observed experimentally in Ref. |43| . the system 
reaches the single-barrier tunneling regime, leading to 
Luttinger liquid behavior in conductance with power-law 
dependence in bias voltage. 

For the sake of convenience, we set the following units 
throughout the rest of the paper: e = h = Do = fcs = 1, 
and the temperature T is in unit of Do = 1. 



B. Useful Mappings 

Our generic model Hamiltonian Eq. ([2J in fact can be 
mapped onto various related model systems as we shall 
discuss below, including the anisotropic Kondo model 
(class (i)), the class (ii) and (iii) systems mentioned 
above. Here, we will address the non-equilibrium trans- 
port through a dissipative resonant level based on one of 



the equivalent models: the two-lead anisotropic Kondo 
model. The mappings for the three classes of models dis- 
cussed below will be derived in an analogous way. The 
general scheme of these mappings is via bosonization fol- 
lowed by re-fermionization (or in the opposite order) 53 ' 54 . 



1. Mapping the dissipative resonant level model 
onto the anisotropic Kondo model 

First, we envision a non-equilibrium mapping reveal- 
ing that the leads are controlled by distinct chemical po- 
tentials. Through similar bosonization and refcrmioniza- 
tion procedures as in equilibrium, our model is mapped 
onto an anisotropic Kondo model^ ' 10 ' 11 with the effec- 
tive (Fermi-liquid) left (L) and right lead (i?) 5 -^ (see Ap- 
pendix A for details): 



H 



K 



fe,7=I/,R,CT=t,4. 



e k ^-y] c ki(7 C kl 



J^s+rlS- + h.c.) 



(4) 



where c][ L , R ^ is the electron operator of the effective 
lead L{R) 1 with a the spin quantum number, 7 = L, R is 
the index for the effective non-interacting fermionic leads, 
5+ = S, S- = d, and S z = Q - 1/2 where Q = dU de- 
scribes the charge occupancy of the level. Additionally, 



'7/3 



y 



a a a8 c k'i36 are the spin- flip operators 
between the effective leads 7 and j3, J^'^ 2 ' cx em- 
body the transverse Kondo couplings, J z cx (1— l/v2o*), 
and /i 7 = ±-jy / l/(2a*). It should be noted that this 
mapping is exact near the phase transition where a — > 1 
or a* — j-^ -> 1/2, and thus /i 7 = ±V/2. Note that 
the above mapping takes a spinless dissipative resonant 
level model with spinless fermionic baths c a —\,2 to the 
anisotropic Kondo model with a "spinful" quantum dot 
(with spin operator given by S + ~ ,z ) and "spinful" con- 
duction electron leads c a 1=L R . The appearance of the 
"pseudo-spin" degrees of freedom in the effective Kondo 
model can be understood in terms of the " charge Kondo" 
effect: the the tunneling between a resonant level (which 
can be represented by a "qubit" or a "spin" ) and the spin- 
polarized leads plays an equivalent role as the "pseudo- 
spin" flips between spin of a quantum dot and that of the 
conduction electrons; and the coupling of the charge of 
the resonant level to the bosonic environment acts as the 
Ising coupling between z-components of the pseudo-spins 
on the dot and in the effective leadsi&. Meanwhile, as 
mentioned above, when the resonant-level model shows 
p-h asymmetry, an additional term tdd^ d appears in the 
Hamiltonian, which is equivalent to a local magnetic field 
acting on the impurity spin tdd^d — > B Z S Z via the iden- 
tification: d) = S+, d = S~, and Sd - 1/2 = S z . For 
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simplicity, we do not intend to study further this p-h 
asymmetry term and focus mainly on the effective Kondo 
model in the absence of magnetic field. Note also that 
the mapping has been derived earlier in Ref. [ID and is 
well-known at equilibrium (Ref. 0). In Appendix A, we 
will provide more details regarding the different theoret- 
ical steps, in particular with a finite bias voltage. 



2. Mapping for a resonant level coupled to a FQHE 

Our analysis for the non-equilibrium transport of a dis- 
sipative resonant level model is applicable for describing 
a resonant level quantum dot coupled to two chiral Lut- 
tinger liquid leads, which is relevant for describing quasi- 
particle tunneling between two Fractional Quantum Hall 
Edge (FQHE) states^ 2 - (the class (ii) model mentioned 
above). In the absence of bias voltage, this case has been 
studied in Refs. IMlil. Via the standard bosonization, 



c Q (0) 



(5) 



the Hamiltonian of such system can be written a o 7 ' 8 ' 11 ? 42 
(see Appendix A.): 



HfQHE = H c hiral + H t + H^, 



(6) 



where the lead term H c hi a i describes two chiral Luttinger 
liquid leads with lead index a = 1,2, Ht denotes the 
tunneling term and the bias voltage term is given 
respectively by: 



H, 



1 



chiral 



2 / „ ^— ' V dx 



H t = t 1 e lipi/VJ? d + t 2 e lip2/ ' /1? d + h.c. 
V 1 



(7) 



where the boson field ip a =\,i denotes the chiral Luttinger 
liquid in lead a, the tunneling between lead and the 
resonant level is given by t a , V is the bias voltage, 
and K refers to the Luttinger parameter. Here, we set 
2na = 1 throughout the paper with a being the lattice 
constant. Through the similar refermionization, we 
arrive at the effective Kondo model as shown in Eq. ( [4J 



with the bare Kondo couplings jj 1 ^ 2 -* 



te 



i(V2- 



)V>2 



J z = 1 — l/y2K. The non-equilibrium RG scaling 
equations for H F q HE have the same form as in Eq. (fT4"f . 



ertheless, the mapping can be straigthforwardly gener- 
alized to the same system but with a spinless quantum 
dot which contains many energy levels. In this case, the 
effective Luttinger liquid parameter K' is modified as: 
±t = y + 1 (see Eq. ( IA22P in Appendix A). More gen- 
erally, the mapping can be further generalized to the sys- 
tem of a many-level (single-level) spinless quantum dot 
with Ohmic dissipation coupled to two chiral Luttinger 
liquid leads with Luttinger parameter K, giving rise to 
the effective Luttinger liquid parameter K defined as (see 
Eq. ( IA24I) in Appendix A): 

K = 4 + 1 + K b (8) 



K h 



(9) 



K 

for a many-level spinless quantum dot and 
1 _ 1 
K ~ K 

for a spinless quantum dot with a single resonant level. 
Details of the mapping is given in Appendix A. 

C. Average current 

We may compute the non-equilibrium current operator 
in the effective models through the mappings. We will 
first compute the current operator within the effective 
anisotropic Kondo model as it is the main focus of this 
paper. From the mapping described in Sec. II. B 1., 
we can establish the invariance of the net charge on the 
resonant level upon the mapping: N\ — N 2 = Nl — Nr, 
where N t = J2ki c \i c ki represents the charge in lead i = 
1, 2, whereas iV 7 = J2k c kia c k'ya- represents the charge in 
the effective lead 7 — L, R. This allows us to check that 
the averaged currents within the Kcldysh formalism 2 ^ are 
the same in the original and in the effective Kondo model 
(see Appendix B for details): 



i[Ql - Qr,H k ] 



)-(H2,I-> R). 



(10) 

Thus, the current / can be computed from the Kondo 
model due to the invariance of the average current upon 
the mapping mentioned above. Note that through the 
various mappings mentioned above, it is straightforward 
to see that the current operator for other related models- 
resonant level coupled to FQHE leads and dissipative res- 
onant level (both small and large in size) coupled to inter- 
acting Luttinger liquid leads-take exactly the same form 
as shown in Eq. (fTUj) . 



3. Mapping for a dissipative resonant level coupled 
to interacting leads 



III. NON-EQUILIBRIUM RG APPROACH 
A. RG equations 



So far, we consider here just a dissipative resonant 
single-level coupled to two non-interacting leads. Nev- 



Now, we employ the non-equilibrium RG approach to 
the effective Kondo model 2 ^ in Eq. In this ap- 



proach, the Anderson's poor-man scaling equations are 
generalized to non-equilibrium RG equations by includ- 
ing the frequency dependence of the Kondo couplings 
and the decoherence due to the steady-state current at 
finite bias voltage^. For the sake of simplicity, we as- 
sume that the resonant level (quantum dot) is symmet- 
rically coupled to the right and to the left lead, t\ = ti 



(or J 



(i) 



r(2) 



J_l). We will discuss in Appendix 



C the more general case with ti ^ £ 2 - The dimension- 
less Kondo couplings as a function of frequency U) exhibit 
an extra symmetry due to the particle-hole symmetry of 
the effective Kondo model: g±( z ){uj) = <?j_( z )(— w) where 
g±( z ){D ) = ./V(0)Jl( 2 ) is the initial value, with ./V(O) 
being the density of states per spin of the conduction 
electrons. Here, we suppress the upper script symbol ^ 
in the Kondo couplings since we will now focus on the 
non equilibrium case V not zero. We obtain^: 



dg z (u) 
d\nD 

dg± (cj) 
dlnD 



/3=-l,l 



(11) 



where 6 W = ®{D— \v + iT\), D < D is the running cut- 
off. Here, T is the decoherence (dephasing) rate at finite 
bias which cuts off the RG fiow-2£. In the Kondo model, 
r corresponds to the relaxation rate due to spin flip pro- 
cesses (which are charge flips in the original model), de- 
fined as the broadening V = T s of the dynamical trans- 
verse spin susceptibility X ( w ) m the effective Kondo 
modeler 



x^M = xo- 



iY s 



(12) 



with X (w) being the time Fourier transform of the 
spin susceptibility x ± {t) = i6{t){[S~ (t), S + (0)]) = 

i9(t)([fi(t)W)jl(0)h(0)]), and X o being x -> = 0). 
Here, we take the pesudo-fcrmion representation of the 
spin operators S + ~' z = Ifa^f 2 7/3 with f a =t,l being 
the pesudo-fermion operator and a + ~' z being the Pauli 
matrices^. 

In the original model the decoherence rate Y corre- 
sponds to the charge flip rates, defined as the broadening 
Yd of the resonant-level (d— electron) Green's function 
(or equivalently the imaginary part of the resonant-level 
self-energy Im(Ed(u))): Y — Yd — im(£d(w)) where the 
self-energy Sd(w) of the d— electron Green's function is 
defined via: 1/G<j(cj) oc (jj + €d + Re(Y>d{u)) + iIm(Y,d{u)) 
with Gd(i-o) being the Fourier transform of the resonant- 
level Green's function G d (t) = i0(t){[d(t), d f (0)]). These 
two definitions for Y agree with each other with the 
proper identification: d = S~ , d^ = S + . 

Note that these RG equations in the present context 
were already discussed in the short Ref. Ha , but now we 
will elaborate the methodology. The configurations of 




FIG. 1: RG flow of g±, cr {V/2) at the transition as a function 
of bandwidth cutoff D (in unit of -Do); the bare couplings are 
g± = — g z — 0.1 (in unit of -Do). We have set V = 0.72 (in 
unit of Do). The decoherence rate F is around 0.00117Do. 



the system out of equilibrium are not true eigenstates, 
but acquire a finite lifetime. The spectral function of 
the fermion on the level is peaked at to = ±V/2, and 
therefore we have g±( z )( U} ) ~ Sx(z)(i^V2) on the right 
hand side of Eq. (jllj) . Other Kondo couplings are not 
generated. From Ref. H3 via the Fermi's golden rule of 
the spin-flip rates Y in the Kondo model, we identify: 

r = J £ Jdu>[n a g 2 z (cj)f u -^(1 - f u _^) (13) 

7,7 / ,er 

+n cr g\{u:)f u - M (l - f u -^,) , 

where f u is the Fermi function. Here, 7 = 7' for the g z {u) 
terms while 7^7' for the g\_{oS) terms with 7, 7' being L 
or R. We have introduced the occupation numbers n a for 
up and down spins satisfying + nj, = 1 and S z = {n^ — 
n\)/2. In the delocalized phase, we get rif = = 1/2, in 
agreement with the quantum Boltzmann equation^!. At 
the KT transition, we can use that g±{uj) — —g z {uj) from 
the symmetry of the Kondo model and that J2a n a = 1. 
Finally in the localized phase, we hav e q \ < — g z , and n a 
satisfies |5 Z | -> 1/2 (see Refs. which remains 
true at a finite bias voltage. 

B. Solutions to RG equations 

Following the scheme of Ref. H3, we solve Eqs. 
(jTTJ) and (fl~3"|) self-consistently. First, we compute 
S-L(z)( w = i^/2) for a given cutoff D. We then 
substitute the solutions back into the RG equations to 
get the general solutions for gx(z)(w) at finite V, and 
finally extract the solutions in the limit D — > 0. When 
the cutoff D is lowered, the RG flows are not cutoff by 
V but they continue to flow for Y < D < V until they 
are stopped for D < Y. 
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FIG. 2: (Color online) F(uj) at T = versus uu across the 
KT transition. The bias voltage is fixed at V = 0.32Do- 
r(w) develops a peak (dip) at uj = in the delocalized (local- 
ized) phase, respectively. At u> « ±V (vertical dotted lines), 
r(w) shows peaks (for localized phase) or dips (for delocal- 
ized phase). Note that T(ld) weakly depends on u for |w| < V , 

r(|w|<7)»r(w = o). 



In FigJT]we show a typical RG flow of g±(V/2) at the 
KT transition as a function of bandwidth D with the 
analytical approximation: g±{V/2) « - ^ s for D > V, 



9±(V/2) 



ior T < D < V, and £_l(V72) 



for L> < T. Here, 2? = Doe 1 ^ 23 ^, with L> being the 
ultraviolet cutoff, and D is the running cutoff scale set 
by the RG scaling equations for gj_/ z . This clearly shows 
that the RG flow of g±(V/2) is stopped at T, a much 
lower energy scale than V. 

Note that the charge (or pseudospin) decoherence 
rate T is a function of frequency, T(oj) in the more 
general and rigorous Functional Renormalization Group 
(FRG) framework 31 . Here, T = T(uj = 0) within FRG. 
Nevertheless, we find r(|o;| < V) at T = depends 
weakly on uj and can be well approximated by its value 
at u> = 0, r(T = 0,w) « r(T = 0,u> = 0) (see Fig.©. 
We have checked that the non-equilibrium current 
/(V, T = 0) and conductance G(V, T = 0) obtained from 
this approximation (r(T = 0,w) « T(T = 0,w = 0)) 
agrees very well with that from the more rigorous FRG 
approach based on the frequency-dependent decoherence 
rate T(T = 0,w) (see Eq. (136]) below) as a consequence 
of the fact that the current and conductance are in- 
tegrated quantities over the frequencies, and they are 
insensitive to weak frequency-dependence of T. In Fig. [3] 
we show the RG flow of the decoherence rate T(lli = 0) 
as a function of D, using the same parameter as in Fig. 
1. One observes that V tends to a finite value as D — » 0. 
The inset shows T as a function of V (see also section 
VI). Note that, unlike the equilibrium RG at finite 
temperatures where RG flows are cutoff by temperature 
T, here in non-equilibrium the RG flows will be cutoff by 
the decoherence rate T, an energy scale typically much 




FIG. 3: (Color online) RG flow of T(D) at T = versus D 
(in unit of Do) at a fixed bias voltage Vo = 0.72Do (vertical 
dashed line) at the KT transition with bare Kondo couplings 
g± = 0.1 = —g z . Under RG, T approaches a constant value 
as D -¥ 0: F(D -¥ 0) w 0.00117D . Inset: F as a function of 
V for the same bare Kondo couplings. 



higher than T, but much lower than V, T < T < V. 
Moreover, T(V) is a non-linear function in V. (For 
example, at the KT transition, T KT {V) oc K/[ln(^)] 2 . 
). The unconventional properties of T(V) lead to a 
non-equilibrium conductance (G{V,T = 0)) distinct 
from that in equilibrium (G(T,V — 0)) near the KT 
transition^. In contrast, the equilibrium RG will lead 
to approximately frequency independent couplings, (or 
"flat" functions g±{uj) rj 5_l, z (^ = 0)). 

Notice that the mapping mentioned above works near 
the KT transition, a* = tt— — > 1/2. However, for a 
general case deep in the localized phase, the effective 



Kondo couplings aquire an additional phase J 



(1),(2) 



ti,2e^^ TwW^ 2 * 1 where the more general form of jj 1 ^ 2 ^ 
and its phase <fi S ;2,i are derived and defined in Eq. (|A10[) 
of Appendix A.. This results in a nonzero bare scal- 
ing dimension^ for jj 1M2) , [jj 1M2) ] = \(^2 - -±=f = 



1- + • This slightly modifies the non-equilibrium 
RG scaling equations to the following form: 



d\nD 

dgA_(uj) 
dh\D 



(3V 



where the linear term 




(14) 



)g± (^) 



Eq. (1141) for g± (u>) comes from the bare scaling dimension 

of jj 1 ^ 2 ^ terms mentioned above, and it vanishes in the 
limit of K — > 1/2, as expected. In fact, this term applies 
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FIG. 4: gx, cr (tj) = —g z ,cr{<jj) at the transition at various bias 
voltages V (in unit of -Do); the bare couplings are g±_ — —g z = 
0.1 (in unit of Do)- The arrows give the values of g±(ijJ ~ 0) 
at these bias voltages. 



to the three models (case (i), (ii) and (iii)) through the 
mappings. Note that the above scaling equations may be 
cast in the same form as in Eq. (1TTT) through redefinition 



of the coupling g z : g z ->• g z = g z 
the previous results remain valid upon the above shift of 
9z- 



IV. NON-EQUILIBRIUM CONDUCTANCE 

In the section, we present our results for non- 
equilibrium conductance. All explicit results will be 
obtained for the KT transition point and the localized 
phase, but not for the delocalized phase. 



Non-equilibrium conductance at the KT 
transition 



At the KT transition, we both numerically and analyti- 
cally solve Eqs. (dJJ) and (in the limit of D -> 0). In 
particular, the approximated analytical solution within 
the approximation Oy ss <d(D — V) due to T <C V is 
obtained: 



p 4 In 

+ Q(V - \w - PV/2\) x 
1 



1 



v 



\u-/3V/2\ 



(15) 



\n[D 2 /V max(|u; - f3V/2\,T)} 41n£ 



The solutions at the transition (denoted g± jC r and g ZtCr ) 
are shown in Fig. SJ Since gj_ cr (u>) decreases under the 
RG scheme, the effect of the decoherence leads to min- 
ima; the couplings are severely suppressed at the points 



I 0.003 - 10 




I(V) via numerics 
I(V)viaEq. (18) 
I(V)viaEq. (17) 
I(V)viaEq. (15) 
Linear fit 




V/D„ 

FIG. 5: Non-equilibrium current at the localized-delocalized 
transition. The naive approximate analytical expression 
Eq. (|lTf) fits well with the numerical result over several 
decades (from V w 10" 12 L>o to V « 10~ 3 L>o, see Inset). How- 
ever, it starts to deviate from the numerical result at higher 
bias voltages. 



u) = ±tt. We also check that g_L, cr (w) = —g z ,cr(^>)- 

From the Keldysh calculation up to second order in the 
tunneling amplitudes, the current reads: 



I=~j duj[j24g±(u}) 2 n a x (16) 
U-^(l-U- m )}-(L^R). 

At T = 0, it simplifies as I = § J^y, 2 dugj_(cj). Then, we 
numerically evaluate the non-equilibrium current. The 
differential conductance is obtained as G(V) — dl/dV. 
The T = results at the KT transition are shown in 
Fig. E]and Fig. © 

First, it is instructive to compare the non-equilibrium 
current at the transition to the (naive) approximation: 



Icr ~ — [ff_L,cr(w = 0)J PS - 



V 



(^(T D /V)) 



(17) 



As shown in Fig. [5l our numerically obtained non- 
equilibrium current fits well with the above analytical 
approximation for V < 0.01D . However, it starts to 
deviate from its numerically obtained values for higher 
bias voltages V > 0.01-Do. This deviation is due to the 
fact that the equilibrium form of the conductance at the 
transition is obtained by treating g cr 3_{uj) a flat function 
within —V/2 < ui < V/2: g cr ±(u) w g cr ±(uj = 0) w 
9cr±(T ^0- We have checked that the equilibrium 
coupling g cr ±(uj = 0) indeed corresponds to g e ^{T = V), 
therefore the transport recovers the expected equilibrium 
form for V — > 0. However, since g±(ui) is not a flat 
function for —V/2 < u> < V/2 (it has two minima at 
u = ±V/2), with increasing V (say for V ~ 0.01-Do) the 
non-equilibrium current exhibits a distinct behavior due 
to the frequency dependence of the coupling. 
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FIG. 6: Non-equilibrium conductance G = dl/dV at the KT 
transition. Go is the equilibrium conductance at the transi- 



tion for T — D : G 
bare couplings g± = - 



g z = 0.1A). 



,T = D ) = 0.005tt with the 



In fact, the more accurate approximate expression for 
the non-equilibrium current at the transition is found to 
be: 
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FIG. 7: Conductance in the localized phase (in units of 
7r). (a) G(V) at low bias follows the equilibrium scaling 
(dashed lines), (b) The conductance G(V)/G C is a func- 
tion of V/T* where we have defined G c = G(a c , V) and 
T* = A>e~' r/ Vsi-9l. 



I(a c ,V) « ^(J[ 5±jCr ( w = 0)] 2 ) 
nV 



(18) 



+ ^ ((1 - f ) [&!,«•(" = ^/2)] 2 ) 



where 



I? 2 



9±,cr(w = V/2) « l/ln(-) 



g±, cr (u = 0) 



(19) 



ln(22? 2 /V r2 ) 41n(D/V)y' 



Here, we have treated g_L !Cr (w) 2 within the interval 
— V/2 < uj < V/2 as a semi-ellipse. 

As demonstrated in Fig. [6j the conductance G(V) ob- 
tained via the approximation in Eq. (1181) fits very well 
with that obtained numerically over the whole range of 
< V < Do. In the low-bias V —> (equilibrium) 

limit, since gj_. cr (w = 0) « j[ e ' cr (T = V) <C 1, we have 

J(a c , (fx.lrC 7 = ^)) 2 ' therefore the scaling of 

G{a c , V) is reminiscent of the equilibrium expression in 

Eq. ©, G(a c ,V) « f (^(T = I/))' = fj^^Tpy. 
This agreement between equilibrium and non-equilibrium 
conductance at low V persists up to a crossover scale 
V ~ 0.01-Drj (determined for the parameters used in Fig. 
[5]). At larger biases, the conductance shows a unique non- 
equilibrium profile; see Eq. (fl"8|). We find an excellent 
agreement of the non-equilibrium conductance obtained 
by three different ways — pure numerics, analytical so- 
lution Eq. (TT51) and the approximation in Eq. (fT5)l . 



For large bias voltages V —> Dq, since 
9±,cr{ UJ ) approaches its bare value g±, the non- 
equilibrium conductance increases rapidly and 
reaches G(a c ,V) ps Go = §.9j_- Note that the 
non-equilibrium conductance is always smaller than 
the equilibrium one, G(a c ,V) < G eq (a c ,T = V), 
since g±(w = ±V/2) < g±(w = 0). Additionally, 
in the delocalized phase for V 3> Tk > 0, the RG 
flow of g± is suppressed by the decoherence rate, and 
Gcx l/\n 2 {V/T K ) (Ref.E3). 



B. Non-equilibrium conductance in the localized 
phase 

In the localized phase, we first solve the equilibrium 
RG equations of the effective Kondo model analytically, 
resulting in 



g£1(t) = 



f( 5 i e ,L(T)) 2 

2cg ± {c+\g z \) 

{c+\g z \?-9l&"D 



(20) 



(21) 



where c = \J g\— g\ . We introduce the energy scale 
T* = Doe^V^-si (which vanishes at the KT transi- 
tion) such that g ( l] loc {T) oc (T/T*f c for T -> 0, leading 

to G«(T) oc (T/T*r. 

At a finite bias, we first solve for the self-consistent 
non-equilibrium RG equations both analytically and nu- 
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FIG. 8: Conductance in the localized phase (in units of 
7r). At large bias voltages V, the non-equilibrium conduc- 
tance G(V) (solid lines) is distinct from the equilibrium form 
(dashed lines). The dot-dashed lines stem from an analytical 
approximation via Eq. (|26[1 , 



merically, resulting in: 
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FIG. 9: Scaling of the conductance G(V) at the KT transi- 
tion (same unit as in Fig. (a). For V>T, the conduc- 
tance follows the non-equilibrium scaling G(a c ,V). (b). For 
V < T, now the conductance follows the equilibrium scaling 
G eq (o!c, T). 



9±,loc (w = V/2) « g ± + ^-[V 2c Vc 2 + A 2 V^ (22) 

2c 

2c 



- ^ A 2 + c 2 ] + — [TVc 2 + B 2 T 2c 



+ 



V c \/c 2 +B 2 V 2c ] 
c 



■In 



C -ln 
2 



BT C + \/c 2 + B 2 Y 2c 



BV C + Vc 2 + B 2 V 2c 

AV 2c + ^ c 2 + A 2yic 



A + ViF+A 2 



g±,ioc (lu = 0)^ 9± + ^-[V 2c Vc 2 + A 2 V^ (23) 
2c 



c 2 
V c \Jc 2 + B 2 V 2c ] 



In 



In 



B ar+ C 2 +B 2 i v ) 2c 



BV C + Vc 2 + B 2 V 2 



AV 2c + Vc 2 +A 2 V 4c ] 



and similarily we get 



[A + Vc 2 + A 2 



A 2 



g z .ioc (lu = V/2)^ 9z + '-^[1~V 4c } 

2c 



(24) 



A 2 

g z ,ioc (w = 0)«j z + -[l 

2c 



V 4 



(25) 



r>2 

— V 2c [l 

c 



-2cl 



where we unambiguously identify A = ^ + ^Qjj . , 
B = AV C ; in this expression, V and T have been 
normalized to Dq. 

The non-equilibrium current in the localized phase 
hociy) is obtained via the same approximation leading 
to Eq. ( ED at the KT transition: 



irV /tt 
~2 
ttV 



(IT 

hoc(V) « — y- \g±,ioc{u 



\- ((! - \) ig±M 



Q)]' 



(26) 



y/2)] 



As shown in Fig. [7J we numerically obtain the non- 
equilibrium conductance in the localized phase. For very 
small bias voltages V — > 0, we find that the conductance 
reduces to the equilibrium scaling: G(V) — > G^ C (T — 
V) oc (V/T*) 4c (see Fig. 0(a) and (b)). For g ±>loc « 
|ffz,ioc| and a* = —> 1/2, we get that the exponent 
4c ~ 2a* — 1, in perfect agreement with that obtained 
in equilibrium at low temperatures: G(T) oc T 2a ~ x 
(Ref. i). At higher bias voltages 0.01D < V < D , 
the conductance now follows a unique non-equilibrium 
form (consult Fig. [8]) whose qualitative behavior is sim- 
ilar to that at the KT transition. Our non-equilibrium 
conductance obtained numerically in this phase is in very 
good agreement with that from the above approximated 
analytical solutions in Eq. [55] (see Fig. [S]). 
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FIG. 10: Scaling of the conductance in the localized phase 
with g± = 0.08, g z = -0.12 (in unit of Do), (a). For V > T, 
the conductance follows the non-equilibrium scaling G(a, V). 
(b). For V <SiT, now the conductance follows the equilibrium 
scaling G eq (a c , T). 

C. Non-equilibrium conductance at finite 
temperatures 



e 




FIG. 11: (Color online) Diagram for the FF current noise 
S(oj). The solid lines represent conduction electron propaga- 
tors; the dashed lines denote the pseudo-fermion propagators. 
The current vertex functions L a p(u}\, CJ2) are denoted by the 
shaded squares. 

been explored experimentally. Therefore, our theoretical 
predictions on the noncquilibrium transport at a dissi- 
pative quantum phase transition offer moltivations for 
further experimental investigations in the regime of our 
interest. The above two qualitatively different behav- 
iors in conductance for V < T and V > T crossover at 
V = T. Note that similar behavior has been predicted 
in a different setup consisting of a magnetic Single Elec- 
tron Transistor (SET) in Ref. [22| where a true quantum 
critical point separates the Kondo screened and the local 
moment phases. 



We have also analyzed the finite temperature profile 
of the non-equilibrium conductance at the transition and 
in the localized phase. We distinguish two different be- 
haviors. At the KT transition, for V > T, the conduc- 
tance G(V, T) exhibits the same non-equilibrium form as 
T = 0, G(V, T = 0) (see Fig. Efa)); while as for V < T 
it saturates at the value for the equilibrium conductance 
(V = 0) at finite temperatures (see Fig. |5fb). In the 
localized phase, while for V < T the conductance satu- 
rates at G(V = 0,T) (Fig. [Ha)), for V > T, however, 
G(V, T) exhibits universal power-law scaling: 

G(V, T)/G{V — 0,T)qc (V/T) 4c (27) 

(see Fig. HuT b)). This universal power-law scaling be- 
havior in G(V, T) looks qualitativly similar to that from 
the recent experiment on the transport through a dis- 
sipative resonant level in Ref. H^. However, these two 
power-law behaviors in conductance at a finite bias and 
temperature are different in their orgins: The authors in 
Ref. |43| studied the quantum critical behavior of a dis- 
sipative resonant level in the regime of the delocalized 
phase (a < a c = 1). As the resonant level is detuned 
from the Fermi level, the system at low temperatures ex- 
hibits power-law scaling in conductance at a large bias 
voltage V > T: G{V/T) cx {^) 2a with < a < 1. 
They showed further that this behavior is equivalent to 
that for a single-barrier tunneling of electrons through a 
Luttinger liquid. By contrast, the Luttinger- liquid-like 
power-law scaling in G(V,T) (see Eq. (|27p) we find here 
is the generic feature of a dissipative resonant level in 
the localized phase (a > a c = 1), which has not yet 



V. NON-EQUILIBRIUM FINITE-FREQUENCY 
CURRENT NOISE 

In addition to non-equilibrium current and conduc- 
tance near the localized-delocalized transition addressed 
above, further insight on the phase transition can be 
obtained from the current fluctuations (or noise). The 
zero frequency shot noise has been used to probe the 
fractional charge of quasiparticle excitations in FQHE 
state tunnelings^ . However, even more useful informa- 
tion can be found in the finite-frequency (FF) current 
noise, which can be used to probe the crossover between 
different quantum statistics of the quasiparticles^S. 
Recently, there has been theoretical studies on the FF 
current noise of a non-equilibrium Kondo dot£i~— . So 
far, these studies have not been extended to the non- 
equilibrium FF current noise of a dissipative quantum 
dot. 



A. Functional RG approach 

To address this issue, we combine recently developed 
Functional Renormalization Group (FRG) approach 
in Refs. 1 1711311 and the real-time FRG approach in 
Ref. |H. Within our FRG approaches, as the system 
moves from the delocalized to the localized phase, we 
find the smearing of the dips in current noise spectrum 
for frequencies ui ~ ±V; more interestingly, we find a 
peak-to-dip crossover in the AC conductance at w « ±V. 
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FIG. 12: (Color online) Diagram for renormalization of the 
current vertex function L a p(ui,Lj[) (the squares). The solid 
lines represent conduction electron propagators; the dashed 
lines denote the pseudo-fermion propagators. Here, the 
Kondo couplings g(io) are denoted by the circles. 

These features are detectable in experiments and can 
serve as alternative signatures (besides conductance) 
of the QPT in the dissipative resonant level quantum dot. 

First, via the above mapping, the current through the 
dissipative resonant level quantum dot is given by the 
transverse component of the current I 1 - (t) in the effective 
anisotropic Kondo model as shown in Eq. (ITOli 16 . Fol- 
lowing the real-time RG approach in Ref. |63j the Keldysh 
current operator through the left lead in the effective 
Kondo model (via a generalization of Eq. (flO)) ) is given 

by: Ii{t): 

x [8+ p (t 1 ,t 2 )Sj(t) + h.c.] (28) 

with a tJ 9 = L,R, S f (t) = f K Ht)af K (t), s^(t u t 2 ) = 
c^(ti)cr ± cj3(t 2 ). Here, L a p{t\—t, t—t 2 ) is the left current 
vertex matrix with bare (initial) matrix elements: L\\ = 
L°rb. = 0, L° LR = -L RL = ig ( [ R = g±, L { [ Z L = L ( RR = g z , 
L°^ R = —L° RL = 0, and k = ±1 being the upper and lower 
Keldysh contour, respectively. The emission component 
of the non-equilibrium FF noise of a Kondo quantum dot, 
S < (t) 1 is given by the current-current correlator: 

S< L (t) = {l£(0)l£(t)) (29) 

Similarly, the absorption part of the noise is defined as: 
S>(t) = (If (t)l£ (0)). Note that the current operator 
lj^{t) is non-local in time under RG; the current vertex 
function L a p{t\ — t,t — t 2 ) therefore acquires the double- 
time structure: it keeps track of not only the times elec- 
trons enter (ti) and leave (t 2 ) the dot, but also the time 
t at which the current is measured 63 . The double-time 
structure of the current operator automatically satisfies 
the current conservation: I^(t) — —I R (t) (Ref. 63). The 
frequency-dependent current noise S(uj) is computed via 
the second-order renormalized perturbation theory (see 
diagram in Fig. [11]). Note that due to the double-time 
structure of the current vertex function L a p(ti,t 2 ), in 
the Fourier (frequency) space, L a p(e + u,e) has a two- 
frequency structure; it depends on the incoming (e + uj) 




FIG. 13: (Color online) 3D plot for Llr(oj\, U12) at zero tem- 
perature in the delocalized phase with bare Kondo couplings 
being g± — 0.05Do, g& = 0.05-Do- The bias voltage is fixed at 
V = 0.32A). 



and outgoing (e) frequencies of the electron (see Fig. [TT|) . 
The result reads: 

S<(u) = J2 -2i?eOD Q/3 H<) (30) 

a,/3=L,R 

where the correlator D a p(w) is computed by the diagram 
in Fig. [UJ 

D a p(cj)< = J U[ X ^(^) X/ (0)]<, 

^-G a (e)G (e + n + Lj) 
x L^(e + u},e)Lj a (e,e + u), 
xm = J ^Gf{u)G f (y + n), (31) 

where G is the Green's function in 2 x 2 Keldysh space, 
and its lesser and greater Green's function are related to 
its retarded, advanced, and Keldysh components by: 

G< = (G K -G R + G A )/2 

G> = (G K + G R — G A )/2 (32) 

The lesser (G < ) and greater (G > ) components of Green's 
function of the conduction electron in the leads and of the 
quantum dot (impurity) are given by: 

G 1/r( 6 ) = iAc&fe-HL/* 
G> /R (e) = iA c (e)(l - f^ L/R ) 

G%{e) = 2mS(e)n frT (e) 

G%{e) = 27ri6(e)(n f(T (e) - 1), (33) 

where A c (e) = 2ttNqQ(Do — e) is the density of states of 
the leads, nfa-(e) = (flfa) is the occupation number of 
the pseudofermion which obeys n/^ + rifi = 1, n/ CT (e — >■ 




FIG. 14: (Color online) 3D plot for Llr(u) 1 ,lj2) at zero tem- 
perature in the localized phase with bare Kondo couplings 
being g']_ — 0.05Do, gi = — O.LDq. The bias voltage is fixed 
at V = 0.32 A, • 



0) = 1/2 in the delocalized phase and n/f(e — > 0) —> 0, 
■^/iC 6 — ^ 0) — ^ 1 in the localized phas o 16 ' 18 . Here, the 
pseudofermion occupation number n/ CT and the occupa- 
tion number on the dot rid are related via (n/-f — n /^) = 
(n d ) - 1/2 (Refs.[Hll3). The renormalized current ver- 
tex function L^p (u>i , oj 2 ) and the Kondo couplings g±(uj), 
g z {u) are obtained from the non-equilibrium Functional 
RG approaches in Ref. HH and Refs. I27ll3ll respectively. 
Carrying out the calculations, the finite-frequency noise 
spectrum reads: 



5<H 



a t jd J 



a,/3=L,R 
X fe— — fe—fifj 



(34) 



where f e -fi a is the Fermi function of the lead a = L/R 



given by / e _ Mc( = 1/(1 + e 
noise spectrum reads: 

1 



). The symmetrized 



S(uj) = ~[S<(uj) + S>(u>)} 



(35) 



with the relation between emission and absorption parts 
of the noise spectrum in frequency space S < (oj) — 
S y (—uj) being used. 

The frequency-dependent Kondo couplings <?j_. z (w) 
and current vertex functions L^p(u>i,uj2) are obtained 
self-consistently within the FRG approaches, which can 
be divided into two parts. First, the Kondo couplings 
g± : z{u) are solved via Eq. ([TT]) 17 ' 27 ' 31 together with the 
generalized frequency-dependent dynamical decoherence 
rate T(ui) appearing in 9 W = Q(D — \co + iT(u;)\) in 
Eq. (jlip . Here, T(uj) is obtained from the imaginary part 
of the pseudofermion self energ y 17 ' 31 ' 57 : 



r( W ) = J Jde g±(e + u,)g ± (e)\f? 

+ g z (e + uj)g z (e)[fL- 
+ (L->R). 



- f R 



L 1 



FIG. 15: (Color online) S(uj) at zero temperature versus uj 
across the KT transition. The bias voltage is fixed at V = 
0.32Do- Inset: S(u) at zero temperature versus w normalized 
to S = S(uj = 0). 



Note that the zero-frequency decoherence rate T(uj = 0) 
corresponds to the decoherence rate T obtained in 
Eq. ([HI 27 . We have solved the RG equations Eq. CCE]) 
subject to Eq. (f3T5)) self-consistcntly^ 5 -. 

The solutions for g±(uj), g z {u) and T(oj) close to the 
KT transition are shown in Refs. [Mlzl (see also Fig. H 
and Fig. @|. As the system goes from the delocalized 
to localized phase, the features in g±(ui) at w = ±V/2 
undergoes a crossover from symmetric two peaks to 
symmetric two dips, while the symmetric two peaks in 
g z (cj = ±V/2) still remain peaks. The finite- frequency 
non-equilibrium decoherence rate T(u>) monotonically in- 
creases with increasing uj, it shows logarithmic sigulari- 
ties at = V in the delocalized phased. As the sys- 
tem moves to the localized phase, the overall magnitude 
of F(w) decreases rapidly and the singular behaviors at 
uj = ±V get smeared oufei 7 .. 

Next, following Ref. HU, we generalize the RG scaling 
equation for the general current vertex function L a p{ui) 
for the anisotropic Kondo model (see diagrams in Fig.[T2l 
and also in Fig. 1 of Ref. l63h . The RG scaling equa- 
tions for the general vertex functions L a 'p Z (u>i,cd2) can 
be simplified as: 



(36) 



dL a p(uji,uj 2 ) , , , , , , 
j/To = L ay {uj 1 ,uj 2 )O ti _ i (uj 2 )g y p(uj2) 

7=L,R 

+ g ai (uj 1 )Q M (uji)L 1 p{uj 1 ,uj 2 ) (37) 

where we make the following identifications: 
9lr/rl{u) -> 9lr/rl{u) = 5J-(w), 9aa(^) 
9ll/rr( u ) = 9z{u). Similarly, L lr/rl (uj 1 ,uj 2 ) -> 
^lr/rl ( w i ' ^2) refers to only the transverse compo- 
nent of the current vertex function L a p(u)\,u)2); while 
Lll/rr ~^ ^ll/rr refers onr y to the longitudinal part 
of Laa ■ Here, the frequency-dependent Kondo couplings 
9±,zo(u) in Eq. (f37"l) are obtained from Eq. (fTTj) and 
Eq. (|3"6")l . Note that the scaling equations for L a p{tJi,u)2) 
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FIG. 16: (Color online) (a) The zero-temperature AC conduc- 
tance Gac(w) defined in Eq. ( I38p versus uj across the KT 
transition, (b) Gac(w) at zero temperature versus uj nor- 
malized to Go = Gac{uj = 0). The bias voltage is fixed at 
V = 0.32A). 



via Ref. |63| can also be expressed within the RG ap- 
proach in Ref. via a straightforward generalization 
by allowing for the two-frequency dependent vertex 
functions L a p(ui\,W2) where U)U2) refers to the incoming 
(outgoing) frequency (see Fig. flTI and Fig. [12]) . 



B. Results 

We solved the self-consistent RG scaling equations 
Eq. (1371) for the current vertex functions with the help of 
the solutions for the renormalized Kondo couplings via 
Eq. pT|) and Eq. (|55|) . The typical results at zero tem- 
perature are shown in Fig. [13] and Fig. 1141 they exhibit 
the following symmetry: L a p(u)i, 0J2) = —Lp a (uj2,uJi). 
Note that since the initial conditions for the current 
vertex function have the following structures: LP aa = 0, 
l lr 7^ 0, we find L aa (uii, ui 2 ) < L lr (loi,lo 2 ). In the 
delocalized (Kondo) phase, a sharp peak is developed 
in Llji(oji,uj2) for (wi,oj 2 ) = (V/2, — V/2); while as 
a small dip is formed for (^1,^2) = (— V/2, V/2). 
Meanwhile, in general Llr(uj\ , UJ2) is maximized at 
uj\(2) = ±V/2 for fixed oj 2 (i)- This agrees perfectly with 
the result in Ref. [fiH In the localized phase, however, 
we find the opposite: Llr(u 1,^2) develops a sharp 
dip at (u)i,u> 2 ) — (V/2,— V/2); and it is minimized 
u>i(2) = ±V/2 for fixed ^2(1)- The peak-dip structure of 
the current vertex function L a p plays a crucial role in 
determining the noise spectrum both in delocalized and 
in the localized phases. 

Substituting the numerical solutions for L a p(wx,w 2 ) 
and g a p(uj) into Eq. (|3~4"]l . we get the zero-temperature 
FF noise S(uS). The results at zero temperature are 
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FIG. 17: (a) £ and (b) F as a function of V/D near the KT 
transition. 

shown in Fig. [15] First, the overall magnitude of S(ui) 
decreases rapidly as the system crosses over from the de- 
localized to the localized phase. This can be understood 
easily as the current decreases rapidly in the crossover, 
leading to a rapid decrease in the magnitude of noise. 
For I a; I > V, S(uj) in both phases increases monotoni- 
cally with increasing uj due to the increase of the photon 
emission at higher energies 6 ^. For \ui\ < V, however, it 
changes from a peak to a dip centered at lo = as the sys- 
tem crosses overs from delocalized to localized phase (see 
Fig. [T"5]) . At |w| = V, S(uj) exhibits a dip (minima) in 
the delocalized phase, a signature of the non-equilibrium 
Kondo effect; while as the system crosses over to the 
localized phase the dips are gradually smeared out and 
they change into a "kink" -like singular point at cj = ±V, 
connecting two curves between u> < V and uj > V . 
We furthermore computed the non-equilibrium AC con- 
ductance at zero temperatur e 63 ' 64 : 



GacM 



S<{u)-S>{u) 



(38) 



across the transition. Note that G(uj = 0) = dl/dV 
corresponds to the non-equilibrium differential con- 
ductance. As shown in Fig. [16] (a), in the delocalized 
phase the splitted peaks in Gac^) at uj = ±V are 
signatures of the Kondo resonant at finite bias, and are 
consistent with the dips at seen in the noise spectrum. 
As the system moves to the localized phase, the overall 
magnitudes of Gac{u) as well as the pronounced 
splitted Kondo peaks at uj = ±V get suppressed; 
they change into dips deep in the localized phase (see 
Fig. \W\ (b)). In response to this change in the splitted 
Kondo peaks, the overall shape of Gac^ — > 0) shows 
a dip-to- hump crossover near uj = 0. Note that the 
suppression of the Kondo peaks for Gac^) at uj = ±V 
corresponds to the smearing of the dips at uj = ±V 
shown in the noise spectrum S(uj) (see Fig. I15p . The 
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above evolution in the noise spectrum matches well 
with the non-equilibrium transport properties studied 
in Refs. ll6lll8L and can serve as alternative signatures of 
the localized-delocalized transition in future experiments. 



To date, the full crossover function of the conductance 
is not known yet. Further study is therefore needed to 
investigate these issues. 



VI. DISCUSSIONS 



VII. CONCLUSIONS 



We would like to make a few remarks before we 
conclude. Firstly, the distinct non-equilibrium scaling 
behavior seen here is in fact closely tied to the non-trivial 
(non-linear) V dependence of the decoherence rate r(V) 
which cuts off the RG flow (see Fig. [T7] (a) and (b)). 
The decoherence rate T near the transition clearly plays 
a very different role as compared to the temperature 
near the transition. In particular, at T = we find that 
r ~ is a highly non-linear function in V, resulting 
in the observed deviation of the non-equilibrium scaling 
from that in equilibrium. In fact, we can obtain the 
analytical form via the approximation in Eq. (1181) and 
Eq. (|26l) . At the KT transition, T/V shows a logarithimic 
decrease as V decreases (see Eqs. (TT9[) . ([20|) ); while in 
the localized phase it exhibits a combined power-law 
and logarithmic dependence on V (see Eqs. ([22)) . ([25]). 

(H, (EH), (EH)). 

By contrast, the equilibrium decoherence rate r(V = 
0,T) shows a clear power-law behavior in the localized 
phase at low temperatures, T ~ > (see Eq. (|21|l): 



T(V = 0,T) cx 



l+4c 



(39) 



which is consistent with the prediction made in Ref. [63 
for the electron lifetime in Luttinger liquids. 

Meanwhile, at the KT transition and in the localized 
phase, since r«V, the RG flow for g±/ z are cut off at 
an energy scale T much lower than V, leading to smaller 
renormalized couplings g±/ z in magnitude compared to 

their corresponding equilibrium values 9±j z (T = V), 

\g±/z\ < \g±j z (T = V)\. This results in smaller conduc- 
tance than that in equilibrium, G(V) < G eq (T). 



Secondly, it is of fundamental importance to study fur- 
ther the possible scaling behaviors in non-equilibrium 
dynamical quantities near the transition, such as the 
lj/T scaling in dynamical charge susceptibility at the 
KT transition and in the localized phase. In particular, 
the question has been raised on the existence of the con- 
cept of "effective temperature" that allows one to extend 
the fluctuation-dissipation theorem in equilibrium to the 
non-equilibrium (non-linear) regime^. It is also interest- 
ing to address the crossover between delocalized phase 
with G(V) cx l/\n 2 (V/T K ) where In TV cx l/(a - a c ) to 
KT point with G(V) oc 1/ In 2 (T/V) and further to the lo- 
calized phase with power-law conductance G(V) cx V? . 



In summary, we have investigated the non-equilibrium 
transport at a QPT using a standard nano-model, the 
dissipative resonant level model. By employing an 
exact mapping onto the anisotropic Kondo model and 
by applying a controlled energy- dependent RG and 
Functional RG approaches to our model system we have 
calculated the renormalized coupling functions g_i_ tZ (u)), 
the decoherence rate F, the current I, differential con- 
ductance G(V,T), and the current noise spectrum S(lu). 
For V — > 0, the conductance G follows the equilibrium 
behavior; by increasing V, the frequency-dependence 
of the couplings begins to play an important role and 
therefore we systematically find scaling behavior of the 
non-equilibrium conductance very distinct from that of 
the equilibrium counterpart. We have also analyzed the 
finite temperature profile of G(V, T) at the transition 
as well as in the localized phase and found that the 
conductance shows different behaviors for V > T and 
V < T; it exhibits V/T scaling behavior for V < T. 

Regarding transport properties of our system near 
the transition, the role played by the bias voltage is 
very different from that played by the temperature. 
The key to these very different behaviors lies in the 
fact that the non-equilibrium charge (or effective spin) 
decoherence rate, which serves as a cutoff for the RG 
flows of the Kondo couplings, is a highly non-linear 
function of the bias voltage. Further investigations 
are needed to address the full crossover function in 
conductance as well as the scaling behaviors of the 
dynamical quantities near the transition in a search for 
the existence of the "effective temperature" that allows 
one to generalize the equilibrium fluctuation-dissipation 
theorem to the non-equilibrium regime. Furthermore, we 
provide signatures of the localized-delocalized transition 
in the finite-frequency current noise spectrum and the 
AC conductance. Our results have a direct experimental 
relevance for dissipative two-level systems; moreover, 
they are applicable for describing non-equilibrium trans- 
port of a resonant level coupled to interacting chiral 
Luttinger liquid generated by fractional quantum Hall 
edge states via the mappings discussed in Appendix A. 
Finally, our model system has direct relevance for the 
recent experiment in a quantum dot coupled to resistive 
environment as shown in Ref. |43T Our work motivates 
future experimental as well as theoretical investigations 
on dissipative quantum phase transitions in nanosystcms. 
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Appendix A: Useful mappings 

In this Appendix, we provide detailed derivations on 
various mappings mentioned in Sec. II. Via bosonization 
and refermionization techniques, the three mappings de- 
scribed below will follow one from the other, but there 
are a few technical details that will change. 



1. Mapping a dissipative resonant level model onto 
anisotropic Kondo model 

We describe in details on the mapping of dissipative 
resonant level model in Eq. ^ onto anisotropic Kondo 
model in Eq. ((4]). Our goal is to connect the parameters 
of these two equations in the main text. 

We first start from Eq. ©: 



# = E ( e ( fc ) _ HiJcliCH + tiC^d + h.c. (Al) 

+ J2 X r( dfd - 1 / 2 )( b r+bt)+J2 UJ rbib r , 



dp 



4>o(x,t) 



) 277^2^1 

fl (x,t) = d t 4> (x,t) 



[b p e lpx + ble* px ]e- aM/2 



(A2) 



where to r — vtp r with Vf, being the phonon velocity, and 
the boson fields cj>o(x,t) and ITo(x,t) satisfies the com- 
mutation relation: [4>o(x, t), IIo(x', t)] = i8(x — x'). The 
dissipative boson bath can therefore be re-expressed as: 



Hdiss = u T blb r = J 7j^\p\blb p 



= - J dx[(d x <f> y(x,t) + Il*(x,t)}. (A3) 

Here, the velocity of the boson field </> is set to be 1. 

We start the mapping by bosonizing the fermionic op- 
erators in the leads: 



2a (0) = 



1 



--F a e 



i<p a (0) 



(A4) 



where we have introduced the (standard) Klein factors 
F a ensuring anti-commutation relations and a is a short- 
distance cutoff (lattice spacing). The fermionic baths of 
conduction electrons can be re-written as: 



Hi, 



E 

fc,i=l,2 



(e(k) - /ii)cLc fe 



= \ [ dx E [(d^ a ) 2 (x,t)+U 2 a (x,t)l^) 



where the Fermi velocity of the electrons is set to be 1. 

The level on the quantum dot can be mapped onto a 
pseudo-spin: d = F d S~ and S z = d)d— 1/2; a = 1,2 
represent the two leads. The coupling between the dot 
and the dissipation bath (A.; term) can be absorbed in the 
tunneling part of the Hamiltonian through the unitary 
transformation Ub^' 



U b 
Hi 



(A6) 



U B H t U B 

J2 toFlFj^^e^S- 



H.c. 



where U is the (real-valued) hopping amplitude between 
the lead i and the quantum dot, Cki and d are electron 
operators for the (Fermi-liquid type) leads and the quan- 
tum dot, respectively, — ±V/2 is the chemical poten- 
tial applied on the lead i (V denotes the bias voltage), 
while the dot level is at zero chemical potential. Here, b r 
are the boson operators of the dissipative bath with an 
Ohmic type spectral density. It proves to be more con- 
venient to re-express the dissipative boson fields b r and 
b\ in terms of the canonical fields (f>o(x,t) and 1X0(2;, t) 



with Kb = —. Here, a refers to the strength of the cou- 
pling between the resonant level and the dissipative bo- 
son bath, and we set 27ra = 1 for simplicity. 

We can simplify our variables even further by combin- 
ing the above fields describing the leads and the noise: 

4>s, a = VK&a + \ Wr4>o), 4>a,a = V^( A 



wher 



be interpreted as the effective Luttinger liquid parame- 
ter as the effect of Ohmic dissipation on the quantum dot 



1 



1 



1 



1 

a* 



Note that here K may 
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plays a similar role as interactions in the Luttinger liquid 
leads coupled to the dot with the identification K = — . 

The comined bosonic and fermionic bath Hbath can be re- 
expressed in terms of these new boson fields: 



+ - J dx[(d x M 2 (^t) + n 2 (x,t)} 

= \ I dx E [(d x 4> s , a ) 2 (x,t) + U 2 Sia (x,t) 

•* a=l,2 

+ (dj a , a ) 2 (x,t) + IL 2 (x,t)}, (A7) 



where H s ( a ) , a fields are canonically conjuate to the fields 
0s(q),q- Note that as we shall see below only the fields 
from symmetric combinations 4> s ,a and Tl s , a couple to the 
tunneling and chemical potential terms, the antisymmet- 
ric combinations <f> a ,a and TL a , a are de-coupled from the 
rest of the Hamiltonian. 

The tunneling and chemical potential parts of the 
Hamiltonian now become: 



H t = U B H t U B 



c*=l,2 



t a FlF d e^S- +H.c(A8) 



+ V / 1 

2 y k 



Close to a* — 1/2 (transition), we can map our 
model onto the 2-channel aniotropic Kondo model. Af- 
ter applying the two unitary transformations Ui — 



e *(^-^.i)^ and u = 



■I v s.Z 
M 7= 



+ t 2 F^F d e iV2 ~VW'^^ e ^^.2]s~ 
- (V2-^=)(dJ s s+dJ sa )S z 



we obtain: 
(A9) 

+ H.c. 



K 

Note that there are additional phase factors 

2 in the hopping terms. Since we are in- 

terested in the physics close to the localized-delocalized 
transition, i.e., K — a* — > 1/2, we may drop these phase 
factors in the following analysis. The chemical potential 
term after the above two transformations now becomes 



Now, we can refermionize the bosons and map our 
transformed Hamiltonian 



Hrlm = Hb a th + H t + 



(All) 



onto the anisotropic Kondo model in Eq. @ via the fol- 
lowing identifications: 

-y/2$.,i = $\-$ R ( Al2 ) 

where F£, R is the Klein factor for the effective lead L 
and R, respectively. 

To see the equivalence between these two models, we 
bosonize Eq. (0} and compare it with Eq. (IA11[) : 



Hk — Hi 



lends 



Hi, 



Hj ± + Hj z , 



J a=L.R 



[(d x <t> a ) 2 (x,t) + n 2 a (x,t)} 



V 

1 



% E ^ 



x^Rl, 



H 



J± 



— J ± LR 



r(2) „+ a- 



J^F^Fie^n-^i, 



H. 



Jz 



= E 



a=L,R 



(A13) 



With the proper redefinitions of the Klein factors: 

F^F^, and the identifica- 



F^ F,, = F n F^ F 1 F 



ntt c4 



tions: d = S~ , d) = S+, d)d - 1/2 = S z , J} Q} = t a , 
J z = 1 we finally establish the equivalence be- 



/2K 



tween a Kondo model with the effective left (L) and right 
lead (i?) in Eq. (0| and a dissipative resonant level model 
in Eq. (jATT]) . 



2. Mapping a dissipative resonant level model onto 
a resonant level coupled to FQHE 



H'l = uluiH^UxUz 



(A10) 



Note that since the hoping H and chemical potential ir„ 

terms involve only SjQ fields, 4> a . a fields are decoupled 
from the Hamiltonian. 



We provide details here on the mapping of a dissipative 
resonant level model Eq. ([2]) onto a resonant level cou- 
pled to Fractional Quantum Hall Edge states (FQHE) as 
shown in Eq. (J5J). 

We start from the Hamiltonian Eq. © describing a 
resonant level coupled to two FQHE states: 



Hfqhe = H c hi 



H t + flu, 



(A14) 
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where the lead term H c hi a i describes two chiral Luttinger 
liquid leads with lead index a — 1,2, H t denotes the 
tunneling term and the bias voltage term H^ is given 
respectively by: 



H 



1 



chiral 



-[-ex:- 



- q=1 ^ 2 \ ax / 



H t = t 1 e^ l/V ^d + t 2 e l ^ / ^d + h.c 
V 1 

(d<pi - d(f 2 ), 



H, 



2 



(A15) 



where the boson field ip a =i,2 denotes the chiral Luttinger 
liquid in lead a, the tunneling between lead and the res- 
onant level is given by t a , V is the bias voltage, and K 
refers to the Luttinger parameter. 

Via similar Unitary transfermations mentioned above, 

m = e *C3r^)S- and U2 = s^~^ s % Eq. © 

now becomes: 

Hfqhe = U\u\HfqheU\Ui = H c } lira i + H t + H^, 

(A16) 

where the tunneling term H t in Eq. (|A15[) becomes (as- 
summing t\ = t% = t): 

H t = t[e l{ ^-^ )v2 e 1 ^ (A17) 
+ e^-^K^e^^S- + h.c. 

- (i-\/^)(5v^i + a^2)5„ 

V ZK 

and the chemical potential term in Eq. (|A15[) becomes 



b. - -I. 



The equivalence between a resonant level coupled to 
FQHE Eq. (j6)) and a dissipative resonant level model 
Eq. © is established by comparing the transformed 
Hamiltonian Hfqhe in Eq. (|A16[) for the former model 
and H RLM (see Eq. ([AT]) . (|ATU)) . and (jATT|) ) for the lat- 
ter one. 



3. Mapping a dissipative resonant level onto a 
dissipative resonant level coupled to chiral Luttinger 
liquid leads 

Below we provide details on the mapping of a large dis- 
sipative resonant level onto a large resonant level (spin- 
less quantum dot) with Ohmic dissipation coupled to two 
chiral Luttinger liquid leads. The mapping is easily ex- 
tended to the latter case with a small (single-level) reso- 
nant level. 

First, we take the same dissipative boson environment 
as shown in Eq. (|A3[) . Via standard bosonization (see 
Eq. ([5])), the Luttinger leads and the chemical poten- 
tial term take the same bosonized form as Eq. (IA5|) and 



Eq. (J7J, respectively. The remaining parts of the Hamil- 
tonian are modified as follows: 

Hdot = H c i + H t + Hdb, 
H d = ^e dk d\d k , 



H t = tJ k ,Jk>S- +h.c, (A19) 

fe,fe',a=l,2 

H db = X r(d{,d k , - l/2)(6r + b\), 
r,k' 

where td k referrs to the energy spectrum of the many- 
levl dot, the electron destruction operator on the dot 
d(Q) is defined as: d(0) = J2kdk, and spin-flip operator 
S represents for the hoping of charge between lead and 
the doti We then bosonize the electron operators in the 
leads (see Eq. ©) and on the dot: d(0) = -±=F d e l<l ' d . 
Via the unitary transformation U b defined in Eq. IA7[ we 
arrive at: 

H t = ulHtU B (A20) 
= £ t a FlF d e i ^°e i ^-^S-+H.c. 

a=l,2 

To further simplify the hoping term, we define new boson 
fields <j> s (a),a y i a linear combinations of the fields (y a (0) 
and 4>d): 

(Pa, a = VK'( -<Pd), 

V K 

(A21) 

with -j£7 — -4 + 1. The combined fermionic baths of the 
leads and the dot are given by: 

Hf = Hi ea ds + Hd 

= \ I dx E l(d^s, a f(x,t)+nl a (x,t) 

a=l,2 

+ (d x <t> a , a f(x,t)+fl 2 a Jx,t)], (A22) 



where H s r a ),a are canonically conjugate boson fields to 
4>s{a),a fields. In terms of the new fields <fis(a),cn the hop- 
ing and chemical potential terms now become: 



H/_i — > Hfj, 

v r~K~ 

= -^\j^[d x {4>as)-d x (cj) aa ) 



(A23) 



We may furthermore combine the boson fields from the 
leads (f> a .a and from the dissipative bath <fta via the fol- 
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lowing definitions: 



Appendix B: Average currents. 



>K{- 



K' 



'K{- 



1 

Kb 
1 



K 1 



(A24) 



where — = + • Upon applying the unitary transfor- 
mation, the combined fermionic and bosnoic baths terms 
become: 

Hf = \j dxJ2 [(d x 4> s ,*) 2 (x,t)+til a (x,t) 

•* a=l,2 

+ (dj a , a ) 2 (x,t)+ni a (x,t)\, (A25) 



where H s ( a ) .a are canonically conjugate boson fields to 



s(a),a 



fields. 



Meanwhile, the corresponding hoping and chemical po- 
tential terms become: 

Ht = t a FlF d eVis- +H.c, (A26) 



a=l,2 



= -^\/¥[^(^i)-^(^,2)]. (A27) 
^ V K 

Via the similar Unitary transformation mentioned above, 



V2<pi)£r„ 



<(- 



Ui = e 2 and U2 = e ^ ' ~ r "'"' " , the tun- 

neling term becomes (assumming t\ — ti = t): 



4=,) 

iJ t = f[e 



= i\/20 s ,i 



(A28) 



+ e ' 'V* )vi e i>&>.s]s- +h.c. 

- (1 - i/X)(0V2^,i + SV20 a , a )5,. 
V 2iv 

The chemical potential term therefore becomes 
Hfj, —> Hp 

= -T vS a *(^.i) " flx(V2^,a)].(A29) 
^ V 2A 

We may now follow the same refermionization procedure 
as shown in Eq. (|A13[) to map our Hamiltonian onto the 
anisotropic Kondo model in the same form as Eq. (j4]) 
with the following identifications: 

jj lM2) oc t a e^-^\ 



J z oc 1 



fi 



1 

2K' 
_ _ V fK 7 
M ~ 2"V X' 



(A30) 



The above mapping can easily be generalized to a small 
quantum dot with single resonant level with K given by 
Eq. @ where the contribution from the many-level big 
dot is absent here. 



In this Appendix, we prove that the average currents 
in the original model I or i is equivalent to that in the 
effective Kondo model I Kondo- The current operators in 
both models are given by: 



d/dt(N 1 - N 2 ) 

iti ^(cj.i rf - d)c k i) - (1 



(Bl) 



2) 



iKondo = d/dt(N L -N R ) 



(B2) 



= iJ^(sl R S+ - s+ L S-) - (1 -+ 2,L^ R) 

On the other hand, from the bosonized forms of the 
two models, (at the transition) we have: 



(Iori) = {d/dt{N x - N 2 )) 



(B3) 



dx{ — [d x cpi - d x tp 2 ]) 



dx \ 2a" ( 1 [ ^ ( ^ m) ~ dx ^ s ^ 



(1 Kondo 

) = (d/dt(N L -N R )) (B4) 
= f dx (d/dt i d ^l - d x $ R }) 

<7=t,l 

= / dx {d/dt[d x (V24> Stl ) - d x (V24> s , 2 )}) 



Therefore, we have 

i^ori) 



'2K 



{I 



Kondo/ 



(B5) 



(or (I orl ) = -jL=(I Kondo )). The above relation ob- 
tained so far from the mapping is exact at finite bias 
voltages. In the limit of our interest K — a* — > 1/2, 

{Iori) {I Kondo)- 

We can also prove this equivalence through Keldysh 
perturbation theory. We now would like to prove that 



where 



{lKondo(t))K — 

6k 



{I Kondo) K — {Iori)t 



1 



(B6) 



(B7) 



x Tr[e-^ HK T C {S^ (-00, oo)I Kondo {t))] 
Z K = Tr[e-P H «T C (S? (-00, 00))] 



S c (—00,00) 



-if c dfHg(t') 



(B8) 
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and 



average current in the Kondo model is given by: 



Z, 

Z„ri = Trie 



-fSH 



T c (S c °"(-°c,oo))] 



5 c °"(-oo,oo) = e^/c^W) 



(B9) 



Here H^(H eq ) is the Kondo (original) Hamiltonian in 
equilibrium (// — 0), T c (- ■ ■ ) orders the operators along 
the Keldysh contour c. 

1. We first show that Zk = Z or i (the two partition 
functions from the original and the effective Kondo mod- 
els are equivalent) where 

Z orl = Tr[e-^ H T c (Sr ; (-oo, oo))] (BIO) 
S c °"(-oo,oo) = e -^Ldt'H^(t') 

To prove this, we first note that the original and 
the effective Kondo models are related by the above- 
mentioned unitary transfromations: Hk — U^HU with 
U = UyfJyUB- The similar relation holds for the cur- 
rent operators: iKondo = U'I or iU. Using the following 
identities: 



oo 

E 

n=0 



/„* cItUHt)H(t)U{t) 
(_l)n rP 



(Bll) 



where 

Tt[A{t)B(t)C{t)\ 



(1t(U\t)H{t)U(t)} 



Tt[C(t)A(t)B{t)] 
Tr[B(r)C(r)A(T)} 



/OO 
J^iG^M ~ <?<») - (1 H- 2, L H- R) 
-oo 



(CI) 



where G£ d (t) = i(s LR Sd(t)). Following Refl58l. the 
Dyson's equation for G^ d (uj) is given by: 



(C2) 



■(•B12) 



- xIr Mx^ M]-(l->2,£->i?) 

where x^ < (t) = (^^(i)), X^O) = 

xj _< = (5j5j"(t)), xi« M, and are time- 

order and anti-timeordered Green's functions, respec- 
tively. The following relations hold among these correla- 
tion functions: 

x <{uj) + X>{uj) = x'M+x'M (C3) 

x > H-x < H = x fl M-x A M 



where X A ( u} ) is the retarded (advanced) Green's func- 
tion, respectively. Straightforward calculation gives: 



X^H = 27r/ w -^(l-/ w _^)5(a;-c(fc)) (C4) 
xJi^M = -27r/ w _ Mn (l -/ u _ M£ ,)(5(w-e(fc)) 



with t = it the imaginary time and A, _B, C being any 
quantum mechanical operators we can then show that 
Zk = Z or i. In other words, when the Hamiltonian is 
under above unitary transformations, the trace in the 
partition function remains unchanged. 
2. In the similar way, we can prove that: 

Tr[e- 0HK T c (Sf (-oo, oo)I Kondo (t))} (B13) 
= Tr[e-? H T c (S c rl (-°o, oo)I ori (t))] 

where we have used Eq. (jBlip and I Kondo — UH or iU . 

From 1. and 2. mentioned above, we conclude 
that (i K ondo(t)}K = (iori(t)}ori holds for all orders in 
Keldysh perturbation theoy. 



The average current reads 



(I) = / [fu-^-U-mfti (C5) 

J — OO 

+ (fl-f a )x<(w) 

where f 12 = 2irpo(Jj^'^) 2 with p being the constant 
density of states of the leads. 

Following Ref. 53, for Ti = Ar 2 , we have 



(I) 



(U-u L - / w -^)f(w)(xfH - X^H) (C6) 



Appendix C: Non-equilibrium current for ii ^ £2- 

In this Appendix, we derive the general expression for 
the average current for t\ ^ From Eq. (|B3I) . the 



where T(uj) = (2npo 



,2 (gjMgjM) 2 



Note that the 



(9iM) 2 +(siM) 2 
Kondo couplings have been genralized to be frequency 

dependent following the noneuilibrium RG approach. 
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